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Abstract. Subsequent to the ideas presented in our previous papers [J.Phys.: 
Condens. Matter 14 (2002) 13777 and Eur. Phys. J. B 42 (2004) 529], we discuss 
here in detail a new analytical approach to calculating the phase-diagram for the 
Anderson localization in arbitrary spatial dimensions. The transition from delocalized 
to localized states is treated as a generalized diffusion which manifests itself in the 
divergence of averages of wavefunctions (correlators) . This divergence is controlled by 
the Lyapunov exponent 7, which is the inverse of the localization length, ^ = I/7. The 
appearance of the generalized diffusion arises due to the instability of a fundamental 
mode corresponding to correlators. The generalized diffusion can be described in 
terms of signal theory, which operates with the concepts of input and output signals 
and the filter function. Delocalized states correspond to bounded output signals, and 
localized states to unbounded output signals, respectively. Transition from bounded 
to unbounded signals is defined uniquely be the filter function H{z). Simplifications in 
the mathematical derivations of the previous papers (averaging over initial conditions) 
are shown to be mathematically rigorous shortcuts. 
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1. Introduction 

Anderson localization |lj has remained a hot topic in the physics of disordered systems 
for a long time (see review articles [3 El II])- Typically the theoretical analysis 
of disordered systems is based either on approximations or on numerical methods. 
However, it is clear that this is cannot be sufficient for systems revealing a phase 
transition e.g. metal-insulator transition as is the case with Anderson localization. 
In such a case an exact solution is greatly needed. 

In our papers [HI [Z| we presented somewhat schematically and by leaving out many 
mathematical details an exact analytic solution of the Anderson localization problem. 
Subsequently there were critical comments |8j that our proof is too short and condensed, 
and more details are needed, the more so, since the mathematical formalism used by us is 
quite new for this scientific community of Anderson localization. This is why we present 
in the present paper the detailed derivation of the analytical solution with illustrations. 
Note that under the exact solution we mean the calculation of the phase diagram for 
the metal- insulator system ^ [7j . We do not calculate transport and other important 
properties. However, the knowledge of the phase diagram permits to understand, how 
these properties can be calculated. It is typical that any exact analytical solution 
is quite lengthy and uses frequently a non-standard mathematical formalism. If this 
presentation is put to the beginning of the paper the physical content of the new theory 
may be easily lost and drowned. A good illustration is the Onsager solution of the Ising 
problem This is why our theory is presented in a series of papers: the mathematical 
idea was presented in [H], with two illustrations: the one dimensional (1-D) case, where 
the solution is well known, and the two dimensional (2-D) case which is non-trivial. 
Several details of the proof were omitted or replaced by simplified arguments, but in all 
cases these simplifications were pointed out. The paper [7j dealt mainly with the results 
for the arbitrary spatial dimension (N-D). 

The structure of the present paper is as follows. In Section 2 we explain how 
the localized states can be treated in terms of a generalized diffusion. This approach 
allows to understand why for defining the phase diagram it is sufficient to solve exactly 
only equations for the joint correlators. We present and solve here the equations for 
arbitrary dimensions. Section 3 deals with the stability of this solution. It is shown 
that the generalized diffusion arises due to the divergence of the fundamental mode. 
The determination of the stability range of this mode permits us to calculate the phase 
diagram. We explain our language of input and output signals, filter function which is 
new for the Anderson localization community. 
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2. Equations for correlators 
2.1. General aspects 

Strictly speaking, the disordered Anderson tight-binding model [T] with Schrodinger 
equation 

J2 tpM' = (E - eM)ipM, (1) 

M' 

where the summation over Ai' runs over the nearest neighbours of site = 
{mi, m2, . . . , m^)}, cannot be solved exactly for arbitrary dimension D. The reason 
is that the random potential Em enters the equation as a product with the random 
amplitude ipM which corresponds to the multiplicative noise case Therefore, the 
exact solution of eq. ((H) discussed in [HI IZ| is possible only under special circumstances 
to be discussed below. 

The fundamental quantity of a disordered system - the localization length ^ - was 
determined in [HI [Zj via the Lyapunov exponent 7. In these calculations the following 
contradictory conditions have to be satisfied. The phase diagram of the system with 
metal-insulator transition should be obtained in the thermodynamical limit |H] (the 
infinite system). However, the determination of the Lyapunov exponent needs the 
introduction of the coordinate system (starting point) which imposes certain limitations 
on the system size. Moreover, a direction for the growth of the divergent quantity should 
be chosen, despite the fact that all space directions are equivalent in eq. ((T]). 

All these conditions are fulfilled for the semi-infinite system, or system with a 
boundary, where the index n = m/) > 0, but all rrij G (— cxd, cxd), j = 1, 2, . . . ,p, with 
p = D — 1. The boundary which is the layer n = defines the preferential direction 
(the axis n) along which the Lyapunov exponent 7 will be calculated. 

It is convenient to interpret the index n = 0, 1, . . . , 00 as the discrete-time, whereas 
all other indices combine in the vector m = {mi, m2, . . . , m.p}. The Schrodinger equation 
can be rewritten as a recursion equation (in terms of the discrete-time, and assuming 
summation over repeated indices) 

V^n+l,m ^n,mV^n,m V^n— l,m ~l~ -^111,111' V^r),,m' ) (2) 

where the operator 

-^^111,111' = E6m,m' — X! '^m",m', (3) 
m" 

is introduced for compactness (summation over m" includes the nearest neighbours of 
the site m). The discrete-time equation Q is a difference equation of the second order, 
which needs two initial conditions. The first natural condition is ipo,in = 0. The second 
initial condition can be presented in the general form ^i^m = without any additional 
conditions for the arbitrary constants a^, except that they are supposed to be finite. 
The rules for the treatment of the field will be explained below, section 13.31 

The recursion equation Q reveals a general feature of the causality - in its formal 
solution ipn+i,m dcpeuds only on the random variables en',m' with n' < n. This gives 
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us a hint for the exact solution of the equation provided that random variables Em 
on different sites do not correlate: {sm^m') ^m,m'- ^^^^ '^^^^ amphtudes ip 
on the r.h.s. of eq. are statistically independent of 6n,m- In other words, while 
performing the mathematical operations in eq. e.g. taking the square of both 

equation sides, with a further averaging over an ensemble of different realizations of 
the random potentials (symbol (...)), the average of the product of amplitudes ip and 
potentials ^ can be replaced by a product of the corresponding average quantities. 
Therefore, the exact solution can be obtained only when the one-sites potentials en,m 
are independently and identically distributed. We assume hereafter existence of the two 
first moments, {en,m) = and (^£n,m) = '^^5 where the parameter a characterizes the 
disorder level. 

2.2. Generalized diffusion 

The ensemble averaging using recursion eq. (jSj) allows us to calculate averages of 
different functions containing the amplitudes ip- The problem arises, which averages we 
can and should calculate? E.g. several even moments of the amplitude were calculated in 
the 1-D case |S| , whereas only the second momentum was calculated for multidimensional 
systems in jHl Ej . It is generally believed that the complete information is available from 
the complete set of all moments jTO]- However, it is not clear, how this information 
should be analyzed. In particular, how can the one parameter of our interest - the 
localization length - be extracted from the infinite set of amplitude momenta? 

It was suggested in |8; to calculate not only the second, and other momenta, 

but also physical values, such as (ln|'?/'|). The average quantities are divided in |E| into 
two categories: containing physical information (e.g. (In and containing no physical 
information (e.g. (I'?/'^)). 

In other words, one has to understand how the choice of a semi-infinite system with 
a selected direction (n axis) and the boundary (layers in the transversal directions with 
n = 0, 1) allow us to detect the localized states under question. 

Let us consider for simplicity a 1-D system, where the amplitude ip = ipn and the 
second initial condition ipi = a. Again, for simplicity let us consider a to be a real 
quantity, i.e. \ip\'^ = ip"^. The recursive equation reads: 



where (e^) = 0, (4) = cr^ 

Let us consider now the random walk problem which is described by the following 
equations 



Considering n as the discrete-time index, both eqs. (jU and (0) describe the 
dynamics of the system with a stochastic time-dependent perturbation 5„. The only 
difference is that the dynamics of the system (jSj) is trivial: ipn = "^o when there is no 
perturbation, 5„ = 0. In contrast, the dynamical system (jD) even for = reveals 



(4) 



(5) 
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proper dynamics (it is a second order equation!). In the band region, \E\ < 2, this 
corresponds to the bounded motion. Therefore, the proper dynamics of both systems 
corresponds to the bounded trajectories, ipn < 
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Figure 1. Power law divergence of second moment for normal diffusion, eq. jSJl, for 
cr = 0.1. 



Assuming now a 7^ 0, eq. (jSJ describes the diffusion. The diffusion motion is 
characterized by divergences, e.g. for the mean time when the system returns to the 
initial state. In particular, the momenta of the amplitude ipn are also divergent with 
the discrete-time n, e.g. 

reveals a power law divergence. Fig. [T] To detect the diffusion, it is sufficient to 
demonstrate the divergence of the second moment of the amplitude and to establish 
its law of time- dependence. Generally speaking, other moments contain additional 
information, which however is not important. The divergence of the second moment 
defines the conditions of the diffusion appearance: if the second moment is divergent, 
so are all other even moments (this is why the choice of a particular moment for 
further analysis is not unique). The law of the time-dependence divergence allows us 
to distinguish between normal and abnormal diffusion. Important is the fact that the 
second moment {ip"^) in eq. © can be found exactly analytically and thus we prefer its 
use (V'^-definition). Speaking formally, the diffusion could also be classified using other 
averaged quantities, e.g. {\4'n\), but such a choice is not convenient for mathematical 
reasons; it hinders an analytical solution. 

As is well-known, in a 1-D system described by eq. (0} any disorder a 7^ leads 
to localization. This manifests itself in the simultaneous divergence, as n —>■ 00, of 
different average quantities of \ipn\, e.g. the linear divergence of (ln|?/'„|) (log-definition 
of localization), whereas the exponential divergence occurs for the powers of [Hj. 
The appearance of the localization in the approach based on eq. (j3)) is equivalent to 
the appearance of diffusion. In eq. © the random perturbation is additive, which 
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Figure 2. Exponential divergence of second moment for generalized diffusion, eq. 
(HI), for a = 0.1 and E ^0. 

determines the linear (power-law) character of the second moment divergence (normal 
diffusion). In contrast, eq. (jH) contains the product of En and ipn (multiplication) 
which determines the exponential character of the divergence. I.e., we can speak of 
the generalized diffusion, analyzing the (ip^) = f{n) divergence as a function of n. The 
exponential localization here corresponds to the exponential growth of f{n) oc exp(27n). 
Fig. 121 where ^ = 1/7 can be interpreted as the localization length {ip"^ - localization 
definition P E]). Note that the ^/'^-definit ion is convenient not only because the 
equations for the second moment can be exactly solved (see below). It is important that 
this definition works also for the non-exponential localization, which corresponds to the 
non-exponential behaviour of f{n). In contrast, the traditional log- definition is valid 
only for the exponential localization, and does not allow an exact analytical solution. 

It was analytically shown ^ for 1-D systems that if the higher moments of the 
random potential e„ can be neglected (for a ^ 0) and we can restrict ourselves to the 
only parameter a, the exponents in all even moments (as well as in log-definition) are 
proportional to each other, they differ only by numerical factors. 

In terms of the generalized diffusion, the results of Molinari p| can be readily 
interpreted as follows. As the disorder increases, all average quantities of {ipnl become 
simultaneously divergent (including powers and Inl'i/'^l). This is indication for the 
appearance of generalized diffusion with the critical disorder cxo = 0. Note that 
in the 1-D case the proper dynamics is unstable and the generalized diffusion arises 
already at infinitesimally small disorder. The proportionality of the different exponents 
(Lyapunov exponents) in the exponential growth indicates that different definitions of 
the generalized diffusion are in fact identical and one can use any of them, e.g the 
■i/'^-definition. 

Comparison of eqs. © and (jH) shows that eq. © describes more complex 
dynamics, where the generalized diffusion can be detected through the second moment. 
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■^n.m)- space dimension D increases, the dynamics turns out to be more stable 

[7j, i.e. the generahzed diffusion arises at a higher disorder level, a > aQ ^ 0. Along with 
the exponential and non-exponential localizations (unstable motion), stable delocalized 
states (stable dynamics) can also arise. The definition of the phase-diagram of the 
system means the determination of the boundaries between the stability and instability 
of the system. 

That is, the two points in [6]: "(i) We are at present not interested in the shape of 
the distribution which is influenced by the higher moments of the on-site potentials but 
in the problem of localization (phase- diagram); (ii) In the analysis of the moments of the 
amplitudes the localization of states finds its expression in the simultaneous divergence of 
the even moments for n oo." mean nothing else but the idea of generalized diffusion. 

Note that the emergence of the generalized diffusion phenomenon and the existence 
of the phase diagrams in dynamical systems with the proper dynamics are analogous 
to critical phenomena in equilibrium and non-equilibrium systems. In this sense, it is 
clear that use of the i/^^-definition for the determination of the phase-diagram is quite 
sufficient. Exactly solvable problems in physics of phase transitions [S], e.g. the Ising 
model, have demonstrated that in spite of the fact that the physical quantities, i.e. long 
range order parameter, magnetic susceptibihty, specific heat, etc. are defined as different 
averages, they still behave similarly, when only the phase-diagram is considered. This 
is not surprising since one deals here with averages obtained for the same statistical 
ensemble. In particular, the critical temperature obtained for the order parameter 
coincides with that obtained for the susceptibility. This occurs for the exact solutions, 
although deviations can arise for the approximate solutions. 

Similarly, the phase diagram for the system with a metal-insulator transition can 
be defined uniquely if one determines any non-trivial average quantity. This is why the 
choice of (In I?/'!) is not convenient. A product of a sum of variables can be represented 
as a polynomial, where each term can be calculated, whereas the logarithm of a sum 
cannot be represented this way and thus the causality principle cannot be used. On the 
other hand, the averages of the powers of the amplitude can be calculated exactly. 

When we speak of the phase-diagram, we mean determination of the boundaries 
of metallic and insulating phases. In this respect, the results [HI Ej demonstrate that 
the second order moments form a closed and linear set of equations and, where the 
only characteristic of the random potentials is the second momentum, o"^, have a simple 
interpretation. The phase boundaries are defined by the only simple parameter a of 
the random potential. That is, calculation of second moments {correlators) is the direct 
way to obtaining the phase diagram. 

2. 3. Correlators 

In general, the field is complex. We can introduce the simplest non-trivial correlators 
based on the ■j/'^-definition: 

x{n)^^i = {ipn^^ipn,i) , (7) 
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y(ri)m,I = l[{rn,m^n^l,l) + (<-l,n.^n,l)] • (8) 

Let us for mathematical simplicity consider real a^- This is also justified by the fact 
that the final eqs. (jHBj) . (jHTf) remain the same. For real definitions can be simplified: 

= {i^n,uii'n,l) , (9) 
y(^)m,l = {lpn,Yn'4'n-l,\) ■ (10) 

Taking square of the both sides of the eq. Q and using the causality principle, one gets 
x(n + = 5m,i(^^x{n)^^^ + x{n - 1)^,1 + (11) 

'^^m,m'a;(^)m',l''^^l',l — '^^m,m'l/('^)m',l — '^^l,l'y(^)l',m- 

Analogously one obtains the additional relation 

y{n + l)n,,i = -yin)i^^ + £m,m'a;(«)mM- (12) 

Note that eq. (fTT|) contains explicitly the diagonal correlators with 1 = m: x{^)m = 
a;('T-)m,m- It is essential that diagonal correlators are positive numbers, x('^)m > 0. 

2.4- Z-transform 

It is convenient to use along the n axis the Z-transform which is common in discrete-time 
systems jTT] : 

X{zU^ = t (13) 

n=0 ^ 

YizU = t (14) 

n=0 ^ 

After transformation one gets 

(z - z~'^)X{z)^^i - a^5m,ix(2;)m " a;(l)m,i = (15) 

'^^m,m'-^(^)m',l''^^l',l " 'Cni,m'^(^)m',l — '^^l,l'^(2;)l',m, 

^^(2:)m,l = -Y{z)l^rn + >^^m,m'A:(z)m',l- (16) 

Hereafter the symmetry properties of the operator are used. 
The solution of the second equation is simple: 

z 1 

.1 ir-^m,m'^(^)m' 1 n -. 

z^ — 1 z^ — 1 

which leads to the equation for X-correlators 



i — — -Cm miX{z)^i I -X{z)^ iiCi/ \, (17) 



2^ + 1 

Z'^ — \ 



m,m' Am]- 

Taking into account the second initial condition, x(l)m,i = cxj^ai. 
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2.5. Fourier transform 

The obtained equation can be easily solved formally via Fourier expansion. Indeed, one 
gets for diagonal correlators 

r fjPk 

xWm = J (^x(^'k)e-'^-, (19) 
x{z,k) = Y.x{z)me''''^. (20) 

m 

A general equation for X-correlators has to be solved using double Fourier expansion: 
X{z, k, k') = E ^(^)m,ie*'''"+'^''. (21) 

m,l 

This leads to 

U{z, k, k')X{z, k, k') = a(k)a(k') + a\{z, k + k'), (22) 

where 

U{z,k,k') = {z- z-') - £^£(k)£(k') + ^^[^'(k) + C\k% (23) 

and 

a{k) = Y.a^e'''^. (24) 

m 

The function 

p 

C{k) = E -2Y,cos{kj) (25) 
i=i 

arises as the Fourier transform of the operator Cm,\- 

2.6. Diagonal correlators 

Assuming k" = k + k', eq. (j22|) can be rewritten as 

U{z, k, k" - k)X{z, k, k" - k) = a(k)a(k" - k) + a'^x{z, k"). (26) 
Using the relation 

x{zX) = J ^X(z,k,k"-k), (27) 

eq.(|26p can be solved elementary for diagonal correlators: 

Xiz,k)=niz,k)x^''\z,k), (28) 

1 ^ o r dPk' 1 



niz,k) ^ J (27r)P f/(z, k', k - k') ■ ^^^^ 

Therefore we have obtained for an arbitrary second initial condition (field am) an 
exact solution for fundamental averages - diagonal correlators. After use of the Fourier 
transform the solution is expanded in modes labelled by the k-vector and describing 
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the oscillations in the solutions in the transversal direction. It is easy to see that 
the ?/^^-definition leads to equations containing only the a parameter of the random 
potentials. If a = 0, all functions 7i(z, k) = 1 and thus x(^) k) = x^^H^j^)- In 
other words, functions x^^H^^ k) correspond to solutions in a completely ordered system. 
Introduction of disorder, a > 0, transforms the initial solution x'^°)(z,k) into x(^;k), 
where the operator T-C{z, k) describes this transformation. 

In our derivation we assumed the system size in a transversal direction to be infinite, 
rrij G {—oo, oo), and thus the obtained equations correspond to the thermodynamic limit 
{L = oo). In the case of finite size of the system, rrij = 0, 1, . . . , L — 1 (cyclic conditions) 
the integrals should be replaced by series. Let us consider for illustration the 2-D case: 
p = D — 1 = 1. Here the index m = m has L values and the number of modes k = A; 
also equals L. That is, analysis of the diagonal correlators leads to L-order matrices 
(vectors). Surprisingly, this obvious fact was interpreted in ref. |Hj that degrees of 
freedom of the initial problem were lost in our paper [H], thus questioning our proof 
of the Anderson theorem. This conclusion was based on the traditional approach of 
the transfer-matrix [HlEl which operates with much larger, x L^, matrices for 2-D. 
The transition from x matrices to L-order matrices was interpreted as an error 
arising due to averaging over initial conditions [H] and thus artificially introducing a 
translation in the transversal directions. However, we do not use in the present study 
such an averaging (see also Section Klfij) . The L-dimensional matrices are the natural 
mathematical tool for the description of diagonal correlators. Moreover, an excessive 
size of X matrices indicates that the transfer-matrix formalism is not adequate. 
This is seen, in particular, from the existence of the so-called trivial eigenvalues of the 
transfer-matrix [HI QHl, which are a- independent. Moreover, the transfer-matrix does 
not permit to make the transition to the thermodynamic limit and thus to calculate the 
phase-diagram (for more details see P2])- These disadvantages are absent in the 1-D 
case, when formally L = and the matrix dimensions coincide. This is why we used in 
jH] the transfer-matrix method only in the 1-D case and exclusively for the illustration. 

3. Anderson localization and stability 

3. 1 . Signals and filters 

The set of equations for the correlators is linear, the same is true for eqs. ()28|) for each 
k-mode. The modes are normal (no mixing of modes with different k values). The 
above mentioned divergence (exponential or non-exponential growth) of the diagonal 
correlators, arising due to localized states, transforms mathematically to the instability 
of set of linear equations. For hnear discrete-time systems under study the most 
adequate formalism is signal theory [TT]. Following [3 El El, let us define x*-°''(^)k) 
as input signals. In our particular case this is a mathematical characteristic of the 
k-mode for an initial ordered systems. (Its sense is explained below). The inverse 
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Z-transform of a given mode can be performed as follows: 

X^'\z,k)^X^^\n,k). (31) 

The input signal is a one dimensional numerical series. The xi^^ k) is the output signal; 
respectively after the inverse Z-transform the latter corresponds to the numerical series 

k), associated with the disordered system, a > 0. The relation between these two 
signals is given by eg. ()28|1 : the output signals are linear transformations of input signals 
performed through 7i(z, k) functions called the system filter ITTj. 

Physical systems with time as independent variable are causal systems [TTj. In 
our problem with discrete-time index n the causality is of primary importance for the 
interpretation of the result. Note that in the general case the inverse Z-transform is 
defined through the complex integral 

h{n,k) = ^ <fn{z,k)z'' — . (32) 
znz J z 

The integration here is performed over the complex plane called the region of convergence 
(ROC) inmi], well-defined for the causal filters. For the causal systems h{n, k) = for 
n < holds, thus after the inverse Z-transform eqs. (j28|) transform into the convolution 
property [TTj: 

n 

xin,k)=Y.hin-l,k)x^'\l,k). (33) 

1=0 

One can see here the linear transformation of input signals into output signals. The 
fundamental property of the signal theory is that the divergence of the output signals is 
related to the asymptotic behaviour of the filter coefficients h{n, k) as n — > oo, which is 
mathematically equivalent to the poles of the 7i{z, k) function of the complex argument 

zm- 



3.2. Fundamental mode 

Since the appearance of localized states leads to the divergence of the diagonal correlator, 
it is necessary to clarify, which fundamental mode k = ko looses its stability and thus 
is responsible for the divergence. Such a problem is quite general in many branches of 
physics. However, in our case it is quite obvious that the fundamental mode is static, 
ko = 0. 

The mode with k 7^ describes transversal oscillations of the diagonal correlators, 
but the sign of the divergent oscillating solution is not defined (see also sect ion IT^ . On 
the other hand, signs of the diagonal correlators are well defined, they are non-negative, 
x(^)m > 0. That is, in our particular case the solution is trivial, the fundamental mode 
is static. It is thus necessary to study the boundaries of the stability of the fundamental 
mode ko = only (the sign of the divergent non-oscillating solution is defined), which 
means uniquely the determination of the phase diagram of the system. 

Assuming x*^°)(2;,0) = S^'^^z), x{z,Q) = S{z), 7i{z, 0) = H{z)), we arrive at 

S{z) = H{z)S^'^\z). (34) 
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After the inverse Z-transform one gets 



Sr. 



Ehn-is^- (35) 

The fundamental input signal S^^^ (z) here and the fundamental filter are defined by the 
integrals: 

^ ^ iz~l)J {2n)P[iz + l)yz-C^ik)y ^^^^ 



H{z) {z-l)J {2n)P [{z + ly/z - C\k)] ' 

Note also the relation 

That is, we have got here a new derivation of relations earlier derived in Refs. [HIIZI- 
Note that for the 1-D case the integral disappears, since p = 0, and one gets 

^^°^(^) = (^^1) [(^+1)2/^-^2]' (39) 

^ = 1 _ ,^2 (^ + 1) 1 .4.1 

H{z) iz-l)[{z + iy/z-E^)y ^ ^ 

3.3. Fundamental input signal 

Let us discuss the physical sense of the input signal signal S^^\z) (or s^^). Note that 
this corresponds to an ideal system with cr = 0. When there is no disorder, = 0, 
the particular solutions of the tight-binding eq. (Q) are bounded functions, plane waves 
exp{ilCA4) with /C = {fci, ^2, • • • , ^d}, provided for a given energy E 

S{}C) = E (41) 

holds, where 

D 

S{iq = Y,2cos{kj), (42) 
i=i 

otherwise particular solutions are not bounded functions, which lie beyond the band and 
have no physical interpretation. If we divide the wave vector into transversal or normal 
directions, /C = {k, kjj}, and keeping in mind that 

/:(k) = 2cos(A;z3), (43) 

for a given transversal mode k the bounded physical solution exists, provided |>C(k)| < 2. 
These relations are sufficient for understanding the physical sense of the input signal. 
Eq. ()36|1 contains in the integrand the function 

(^ + 1) 1 



{z-l) {z + lY/z-C\\i)' 



(44) 
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Under the condition |>C(k)| < 2 its inverse Z-transform gives 



(45) 



sin (ko) 

where k^, = k£,(k) is the solution of eq. pHjl . Assuming that the arbitrary field a^n is 
selected in such a way that the Fourier transform coefficients a(k) are not zero only if 
|£(k)| < 2, eq.(jnni) after the Z-transform reads 

(0) ^ r dPk \a(k)\^ sin^kpn) 
J {27t)p sin^ (kD) ' 

The s^"^ is bounded for any n. 

On the other hand, assuming that |>C(k)| < 2 is violated for certain k, so the 
a(k) 7^ for |>C(k)| > 2, the inverse Z-transform of eq. ()44|) leads to the function 

sinh^iKpn) ^^^^ 

which increases without bounds as n — >■ cxd. Here Kp = K,D(k) is the solution of equation 

|£(k)| = 2cosh(/tz)). (48) 

The input signal s^^ is also increasing to infinity as a function of n. Note that both 
input and output signals are real values. 

In other words, on can conclude that any physical solution of an ideal system (cr = 0) 
with bounded wave functions has a one-to-one correspondence to the mathematical 
object - a bounded 1-D input signal Such solutions always exist, if the energy 

E lies within the band interval. On the other hand, formal solutions without physical 
interpretation {unbounded solutions beyond the band) correspond to the 1-D unbounded 
input signals. Note that particular numerical values of input or output signals are not 
important from the point of view of the signal theory [2] whose main concern is a 
qualitative discrimination between bounded and unbounded signals. 



3.4- Fundamental filter 

The quantity characterizing the phase diagram of a disordered system is the fundamental 
filter, eq. ljHTj) . Its idea is quite clear Let us consider the physical states inside the 

band, \E\ < 2D. All these states are described by wave functions bounded in amplitude, 
which corresponds to a bounded 1-D input signal s^^-* . In a disordered system this signal 
transforms into s„, according to eq. (jHSj) . This output signal s„ can be either bounded, 
or infinitely growing (generalized diffusion). 

The natural interpretation of bounded output signals is that they - as before - 
correspond to the delocalized states bounded in amplitude. The unbounded output signals 
correspond to the localized states, respectively, which lead in the semi-infinite system to 
the divergence of the diagonal correlators. An important result of signal theory JT] is 
that the divergence does not result from properties of a particular bounded input signal; 
the cause lies in the filter, H{z) or hn- We define the phase diagram of the Anderson 
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model based on a general concept of the signal theory known as a BIBO stability |llj . 
Namely, a system is BIBO stable if every Bounded Input leads to a Bounded Output. 
A stable system (delocalized states) is characterized by a stable filter H{z); its main 
property is the absence of poles in the complex z-plane outside the circle \z\ > 1. An 
unstable system (localized states) is characterized by the unstable filter H{z) with poles 
outside the circle 1^1 > 1. As was shown in [HllZl, in the Anderson model the pole zq 
of the unstable filter lies on the real axis. Consequently, zq = exp(27), where 7 is the 
Lyapunov exponent (■i/'^-definition) |S1 E! • Using the inverse Z-transform, one can easily 
find that asymptotically /i„ ~ exp(27n), i.e. it is divergent. The localization length is 
defined as ^ = 7"^. 

The mathematical formalism of the signal theory allows to extend and complement 
this result for the energy range \E\ > 2D, where in the absence of disorder there existed 
only mathematical (divergent) solutions having no physical interpretation. It is known 
that disorder extends the band, i.e. new physical states arise also at \E\ > 2D. The 
question is: can delocalized states exist amongst these new states? The answer is simple. 
In the region \E\ > 2D only unbounded input signals s^^-* exist. Such a signal cannot be 
transformed into a bounded output signal s„. The output signal s„ is always unbounded 
with a dual interpretation: on the one hand it corresponds to mathematical solutions 
(no physical interpretation); on the other hand, the divergence of the output signal can 
be associated with an emergence of the localized states. I.e., an emergence of delocalized 
states in the region > 2D (outside the band) is impossible [7j. 

Use of the filter function H{z) is a general and abstract method for describing the 
metal-insulator transition, valid for any space dimension. Earlier [HI we used a less 
fundamental approach for the 1-D problem, based on the calculation of the Lyapunov 
exponent for the ?/^^-definition 

7 = lim 7^1ns„, (49) 

where for the 1-D case s„ = (i'n)- Taking into account that the divergence of the output 
signal Sn is caused by the divergence of the filter hn, one can use also the following 
equation 

7 = lim ^Inhn. (50) 

In the 1-D case both definitions of 7 are equivalent, since for a fixed energy and disorder, 
for a second initial condition ipi = a, one gets the single output signal s„. 

Let us use for illustration eqs.(|H^. (jlUj) . Restricting ourselves by the band center 
E = (which simplifies the equations) and making the inverse Z-transform, one gets 

Sn = \a\'^^^^[zi;-{-l/z,)% (51) 
hn = Sn,o + rlxTTt^o - (-1/^0)"], (52) 
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provided s^^ = |ap[l — (— 1)"]/2 (bounded input signal). Indeed, simultaneous 
divergence of the output signal s„ and the filter hn arises due to the fact that 
the asymptotic behaviour of both quantities is determined by the same parameter, 
Zo = exp(27) > 1. 

However, this equivalence is no longer valid for a space dimension higher than 
one. Since the second initial condition is defined by the field a^n, this corresponds 
to a continuum of input signals s^^ and, respectively, a continuum of output signals 
Sn- That is, the definition is not valid, since it is not clear that different signals 
should correspond to the same Lyapunov exponent 7. Formally, a whole continuum 
of solutions should be analyzed. However, eq. ()H5j] demonstrates that the fundamental 
Lyapunov exponent 7 does not depend at all on the field a^, it is sufficient to define 
the Lyapunov exponent using only eq. ()5n|l . 

3.5. Phase diagram and multiplicity of solutions 

A more detailed analysis 01111121 reveals another problem of the definition (jl^ : this 
is valid only in the 1-D case where for any disorder a only the phase of the localized 
solutions exists. 

As it was mentioned above, for the calculation of the fundamental filter hn using 
the inverse Z-transform, eq. ()H8|1 . the contour integration over the so-called region of 
convergence (ROC) JJJ is necessary, provided the filter under consideration is causal. 
It is shown ^ ^ that dependent on the energy E and disorder a, two general cases are 
possible for the Anderson model with D > 2. In the first case, the ROC lies outside 
the circle \z\ > 1. The filter hn is thus unstable and describes the localized states 
(insulating phase). In the second case, the ROC consists of two domains in the complex 
plane: one domain is inside the circle \z\ < 1, another one - outside this circle, \z\ > 1. 
Respectively, there are two ways to calculate the integral using eq. The double 

solution arises, one describing delocalized states (filter hl^\ metallic phase), the other 
localized states (filter hl^\ insulating phase). 

Consequently, in this range of parameters E and a the two phases can co-exist, this 
is why the metal-insulator transition in the Anderson model has to be analyzed in terms 
of first- order phase transition theory IHIIZIE]- For arbitrary random potentials the wave 
function can be either localized or delocalized (no co-existence!). However, in terms of 
statistics of an ensemble of random potentials, both localized and delocalized solutions 
can arise with a comparable probability. Namely this comparability of probabilities is the 
main characteristic of the first-order phase transition, in contrast to the second-order 
transitions where either a pure metallic phase (no localized states) or a pure insulating 
phase (no delocalized states) should exist: co-existence is impossible. Strictly speaking, 
the type of phase transition in the microscopical Anderson model cannot be determined 
without exact solution of the problem. Despite the fact that first-order phase transitions 
are very common in nature, its realization in the Anderson model (under strong influence 
of the phenomenological scaling theory of localization ^21) never was seriously discussed. 
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Phase co-existence creates a serious problem of the choice of an adequate 
mathematical formalism. Strictly speaking, when calculating the average quantities on 
the ensemble of different realizations of random potentials, the contributions from pure 
metallic and insulating phases are considered as equivalent. However, such heterophase 
averages have no physical sense [HI El ^] ■ In a two-phase system with first-order phase 
transition one is interested in properties of pure phases. From this point of view, the 
earlier introduced output signals s„ are also heterophase averages. This is why the 
derivation of equations for signals is not a goal, but the tool for obtaining a more 
fundamental property - the filter H{z) {hn). As was noted [3 El (121, it is such a filter 
which reveals the multiplicity of solutions and which permits us to determine the phase 
diagram for the Anderson model. 

3.6. Averaging over initial conditions 

Let us discuss now the problem of averaging over intitial conditions [HI EI. This 
procedure was introduced in Ref. its meaning is quite simple. First of all, 

we choose a particular second initial condition, i.e. a field a^n- Secondly, let us 
consider the field a'^ = am+mo obtained from the first field as the result of a trivial 
translation in transversal direction by the vector mo. It is obvious that the relevant 
diagonal correlators x'(^)m = x(^)m+mo can also be obtained by the argument shift 
in a transversal direction. Both solutions are physically equivalent. The Fourier 
transform of the field a^i with the vector shift mo satisfies a simple relation: a'(k) = 
exp(— ■ikmo)Q;(k). 

Eqs. (j^Hj) to (jHUj) clearly demonstrate that the vector shift affects only signals for 
particular modes k, and all physically equivalent solutions differ from each other only 
by phases: 



Note that only the fundamental mode k = remains invariant. Averaging now the signal 
over all translations in transversal directions (averaging over initial conditions [Hj) gives 
zero for all non-fundamental modes, k 7^ 0: x^'^K^^'^) — 0) k) = 0. Nonzero is 
only the fundamental mode (k = 0): x^^K^^^) = S^^\z), xl-^, 0) = S{z). The average 
diagonal correlator loses its dependence on the argument m: x("-)m = That is, 
averaging over initial conditions is an efficient tool for getting rid of all non-fundamental 
modes which are non-essential for the analysis. 

Note that eq. (jHEj) contains r(k) = |a(k)p which is a function of the field 
am, its Fourier transform is T^. The averaging over initial conditions corresponds 
mathematically in the initial eq. (fTSj) to the replacement of the initial condition 
x{l)m,\ = ttmtti by the average, x(l)ni,i = Fm-i. After such a replacement of initial 
conditions the system becomes much simpler: it is translation-invariant in transversal 
directions which makes use of the double Fourier transform unnecessary. In other 
words, averaging over initial conditions is nothing but a simple mathematical trick which 




(54) 
(55) 
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permits to get quickly the fundamental property of system - the filter H{z). Such tricks 
are based on the fact that the filters are universal system characteristics, describing 
the solution transformation from order to disorder. Such universal characteristics 
are independent of the initial conditions and other details. The initial conditions 
determine such non-universal properties as signals. Consequently, in order to obtain the 
fundamental filter H{z) and then the phase diagram, one can perform linear operations 
with signals, in particular, averaging over initial conditions. 

Another example of such operations is the following. Let us consider the field ctm 
not as fixed (second initial condition), but as a random variable, with simultaneous 
averaging over the field am and over the ensemble of the random potential realizations 
Em- We assume the absence of correlations between the potentials and the 
field a^n] unlike the correlation between the field components characterized by the 
arbitrary correlation function (am^i) = Tm-i- The initial condition in eq. ()18j) becomes 
x{l)m,\ = Tm-i- Taking into account that the operator £m,m' by its definition, eq.©, 
depends only on the on the argument difference, m — m', one gets the translation 
invariant system, where = X{z),J^_l, provided for the diagonal correlators 

X{z)m,m = x{^)m = X{z)o = S{z) holds (the diagonal correlators are m-independent). 

In this case in order to solve eq. (fTHj) , it is sufficient to use a single Fourier transform 

X{zU = J ^X(^,k)e-'^-. (56) 

Instead of eq.(|^ one gets 

U{z, k, -k)X{z, k) = r(k) + a^S{z), (57) 

which gives 

^(-) = /(^^(-'k). (58) 

As a result, one returns to the fundamental eqs.dHH) to (jTTj) . with replacement of the 
|a(k)|2 for r(k). 

4. Conclusion 

We would like to stress that in this paper we presented a mathematically rigorous method 
for the calculation of the phase diagram for the Anderson localization in arbitrary 
dimensions, which was briefly discussed earlier [HI U\- The phase diagram for the 
metal-insulator transition is obtained using the Lyapunov exponent 7. Localized states 
correspond to values of 7 > 0, i.e. a divergence of the averages over wavefunctions. 
This divergence is mathematically similar to the divergence of averages for the diffusion 
motion. That is, transition to the localized states can be treated as a generalized 
diffusion. From this viewpoint, in order to determine the range of the existence of 
localized states (i.e. the phase diagram) and the type of localization (exponential or 
non-exponential), it is sufficient to solve equations for the joint correlators. We have 
shown that these equations are exactly solvable analytically. 
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In its turn, the appearance of the generahzed diffusion arises due to the instabihty 
of a fundamental mode corresponding to correlators. The generalized diffusion can be 
described in terms of signal theory, which operates with the concepts of input and output 
signals and the filter function. Delocalized states correspond to bounded output signals, 
and localized states to unbounded output signals, respectively. Transition from bounded 
to unbounded signals is defined uniquely by the filter function H{z), or more precisely, 
by the position of its poles in the complex plane. This function can be calculated for 
arbitrary space dimension D. 
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